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Abstract 

Using Sphinx tt, a new MC simulation program for transverse 
polarized nucleon-nucleon scattering based on Pythia 5.6, we cal- 
culate the transverse double spin asymmetry A TT in the Drell-Yan 
process. If one assumes (quite arbitrarily) that the transversity par- 
ton distribution 5q(x,Q 2 ) equals the helicity distribution Aq(x,Q 2 ) 
at some low Qq scale, the resulting asymmetry is of order 1%. In this 
case is A TT would hardly be be measurable with PHENIX at RHIC. 
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1 Introduction 



Originally, the transversity parton distribution 5q(x, Q 2 )[]was defined by Ral- 
ston and Soper Q. However, only after polarized targets were developped 
and first spin-experiments yielded very interesting results, transverse spin 
structure functions became a real issue. After the "rediscovery" of Sq(x, Q 2 ) 
by Artru and Mekhfi || in 1990 the topic of transverse spin has gained an 
increasing amount of attention by elementary particle physicists. Much of 
the interest is caused by the fact that besides the well known unpolarized and 
helicity weighted parton distributions q(x, Q 2 ) and Aq(x, Q 2 ) the transver- 
sity distribution 5q(x, Q 2 ) is the only twist- 2 distribution of the nucleon. On 
the other hand, its special chirality selection rules prevent its measurement 
at leading order in polarized inclusive deep inelastic lepton-proton scattering 
experiments and therefore no experimental data on Sq(x, Q 2 ) is available up 
to now. This situation is expected to change in the near future, when the 
Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory 
will be able to scatter transversely polarized protons with high cm energy, 
luminosity, and large polarization of each beam. One of the proposed ex- 
periments for the RHIC spin program is the measurement of the transverse 
double spin asymmetry 

daTT + daTi {) 

using the Drell-Yan process q^ — > 0, |, |5[], where i denotes the flavour 
and an electron-positron or muon pair. Here, the asymmetry is not 
suppressed because the two quark lines in the unitarity diagram are inde- 
pendent of each other 0. This also holds for other partonic processes, like 
the prompt photon production with an away-side jet qi + qi^ jg, or two-jet 
production qi + q\ — > gg, qiqi — > qiqi. Nonetheless, these reactions are not 
very useful for extracting information about the transversity distribution, 
since either the partonic part of the asymmetry is small or the cross section 
of other contributing background processes is much larger || |7j . 

It is generally believed, that at some low scale Qq the transversity 
distribution is similar to the helicity distribution 

5q(x,Q 2 )^Aq(x,Q 2 ) (2) 



1 Othcr common denotions for Sq(x,Q 2 ) are hl(x,Q 2 ), h^(x,Q 2 ), and Aiq(x,q 2 ). 
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and that their difference is only due to dynamic effects || |9|] . Crude analytical 
estimates which were based on this assumption yielded A TT ~ 10% || |Hj] . 
However, in these model calculations the correct Q 2 -evolution of Sq(x, Q 2 ) 
was not taken into account. 

In this paper we present a new estimate of the transverse double spin 
asymmetry of the Drell-Yan process based on (0) using the new MC simu- 
lation program Sphinx tt. We take into account the correct Q 2 -evolution 
and give a crude estimate for the QCD corrections to the simple Drell-Yan 
model. 

In section ^ we describe the model we used for our calculations. Sec- 
tion || contains a short description of Sphinx tt. Our results and conclusions 
are presented in section 

2 Our model 

In order to be able to make an estimate of A TT we assume that the 
transversity weighted distribution of quarks and antiquarks equals the he- 
licity weighted distribution at the low scale Qq = 4 GeV 2 : 

5q(x, 4 GeV 2 ) = Aq(x, 4 GeV 2 ), 8q(x, 4 GeV 2 ) = Aq(x, 4 GeV 2 ) . (3) 

This is rather an ad hoc assumption. Transversity distributions and helic- 
ity distributions should agree for some very small Q 2 , where perturbative 
QCD is not applicable. It is not clear how much they will differ after evo- 
lution up to 4 GeV 2 . Their DGLAP equations differ most notably because 
the transversity distribution is not a singlett with respect to chirality and 
therefore does not couple to gluons. Its evolution equation reads simply 

_^ (I ,o. ) = ^) J f;^ (y , < n,p fI (£) . (4) 

The splitting function 8P qq {z) = C 2 {R)[2/(1 -z) + -2 + (3/2)5(z - 1)] was 
calculated by Artru and Mekhfi 0. Equation (|j) implies for the anomalous 
dimensions 

H W (Q 2 ) = f'dxx^Sq^Q 2 ) , (5) 
Jo 
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SP$f> = l dzz n - 1 6P qq (z) = 0,^)^-2^- J , (6) 

( n {r a\\ SPg>/2xb 
^ (n) (Q 2 ) = Sqt\Ql) MM) . (7) 



Since b = (33 - 2n f )/12n > 0, < and 5P$+V < &p0O all mo- 

menta of the transversity distribution decrease and the distribution becomes 
softer with increasing Q 2 . 

We used the leading-order parametrizations Set C of Gehrmann/Stir- 
ling |12[], and Bartelski/Tatur |L3| and evolved them according to 



The authors of all major previous works fl], §, |], §, [1]] about the trans- 
verse double spin asymmetry of lepton pair production only considered the 
simple electromagnetic graph qq — > which yields the hadronic double 
spin asymmetry 

TT = sin 2 6 cos 20 gj e?5qi(x A , Q 2 ) Sqj(x B} Q 2 ) + A <-» B 

DY ~ 1 + C0S 2§ Ei ehi( x A, Q 2 ) Qi(xB, Q 2 )+A<r->B ■ U 

9 and 4> are the polar and azimuthal angles of the lepton momentum in the 
dilepton rest frame with respect to the beam and spin axis, resp. 

However, it is well known that the cross section of the QCD corrections 
at 0(a s ) is approximately as large as the one of the Drell-Yan graph resulting 
in a f^-factor of 

k dy = dg o(a 0) + da 0(as) ^ 2 ^ 
doei(aO) 

Since this Bremsstrahlung graphs might carry asymmetries completely dif- 
ferent to @ one should have a closer look at them. In general, one would 
expect the analytical expressions for these asymmetries to be more compli- 
cated than (||) since one more outgoing parton is involved here. So it would 
be nice, if one could seperate the Drell-Yan graph from its QCD corrections 
kinematically. 

The transverse momentum of the virtual photon produced in the 
Drell-Yan model is of the order of 1 GeV since the quark and antiquark 
possess a primordial transverse momentum in the proton of the order of 
several hundred MeV. Contrarily, the virtual photon in the QCD-Compton 
and QCD-annihilation process can also posses a large transverse momentum, 
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(a) (b) (c) 

Figure 1: QCD corrections at order a s to the simple Drell-Yan model: (a) 
annihilation graph (b) Compton graph, (c) virtual corrections. The decay of 
the virtual photon 7* is not shown. 



which is balanced by the large opposite transverse momentum of the ad- 
ditional outgoing parton. In case the transverse momentum of the virtual 
photon in the QCD annihilation or Compton process is small, the reaction 
is characterized by three scales, namely q 2 , q^, Aq CD . Here, the differential 
cross section is dominated by double-leading-logarithmic contributions 

b l + b 2 a s {q 2 )\n 2 l^i+ha 2 s {q 

9 ^ ' (10) 

where the hi are constants. The true expansion parameter is rather 
a s (q 2 ) \n 2 (q^/q 2 ) than a s {q 2 ). Equation (|T0|) describes the radiation of an 
infinite number of soft gluons from the partons which enter the hard par- 
tonic reaction. Fortunately the coefficients hi are not independent of each 
other but exponentiate |L4] . 





The differential cross section da/dg^dg 2 at double-logarithmic precision 
differs from the one calculated by a naive convolution only by a form factor Tp 
for the Q CD-annihilation process and T f Tq for the QCD-Compton process 
where |Tj| 

T F (q 2 ,q 2 ) = e- B / 3 , T G (q 2 ,q 2 )=e-^\ B = (a./*) In 2 {^\ . (11) 

These form factors make (render) the differential cross sections finite at qr = 
0. According to EB1 they also contain the virtual corrections. 
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If we assume, that the partons inside the proton contain an intrin- 
sic transverse momentum which is distributed according to a Gaussian 
dN/dk T oc expk^/a 2 with (k T ) ~ 600 Me V then the theoretical predic- 
tions describe the data well |14[]. This value is perhaps somewhat high, as 
deep inelastic lepton-nucleon scattering experiments suggest a smaller value 
of about (k T ) ~ 400 MeV. 

Actually, the double-logarithmic approximation (DLA) does not repre- 
sent the "state of the art". If one incorporates sub-leading terms [[17]] then 
primordial momentum is not needed to describe the data. This implies that 
the contribution from the Drell-Yan graph vanishes because da/dg-r = at 
qr = and the resulting shape of the 0(a s ) contributions is rather changed. 
Such an approach would go beyond the present calculations and would be 
required if our calculations suggested a statistical accuracy to be expected 
from the planned RHIC experiments (which they do not). Thus, we restrict 
ourselves to the DLA for the calculation of the event rates although strictly 



speaking it is not valid for all values of qr |T7 



After having clarified the situation for the unpolarized cross section 
we want to consider polarized scattering now. Transverse spin asymmetries 
manifest themselves in a nontrivial distribution of the azimuth of the out- 
going leptons (see (||)). Their origin can be easily understood. The axis of 
the incoming momenta in the cm frame and the spin axis define one plane 
in space. The orientation of the second one, defined by the momenta of 
the outgoing leptons and the ones of the entering partons, is therefore not 
arbitrary. 

Figure |]a shows the situation of the unpolarized QCD-annihilation 
graph which is similar to the previous one. The orientation of the second 
plane in relation to that of the first one is not arbitrary again. The major 
difference to the previous case is that the first plane, defined by the momenta 
of the incoming partons and the gluon, now is not fixed in space because the 
azimuth of the gluon momentum is not fixed. The same situation also holds 
for the QCD-Compton process. 

Unfortunately, the calculation of the resulting angular distributions in 
QCD is not well understood. One can try to use only perturbation theory 
although qx becomes very small in the region we are interested in, or one can 
invoke resummation techniques for the small qx region but neither approach 
describes the data completely JT8| . 

In case the incoming quark and antiquark are transversely polarized, 



6 



an additional third plane can be defined according to figure 0b. Here, the 
orientations of all planes in relation to each other are nontrivial. The virtual 
photon and the gluon now take on similar roles as the lepton pair in the Drell- 
Yan graph since they are the outgoing particles from the first interaction. 
This means, that the distribution of the azimuthal angle between the photon 
momentum and the spin axis shows the typical cos 2(f) modulation. 

If one of the entering particles is a gluon then the azimuthal angle be- 
tween the first and the third plane is distributed evenly in the allowed inter- 
val, since gluons cannot carry an asymmetry. For this reason, the Compton 
graph does not carry a transverse spin asymmetry at all. 

We did not make an attempt to calculate the transverse spin asymmetry 
of the QCD-annihilation process because, as we mentioned above, the angular 
distributions in the unpolarized case are not well understood. Thus, one 
might expect to run into similar difficulties in the polarized case. 




(a) (b) 

Figure 2: (a) the unpolarized QCD-annihilation process, (b) the case with 
transverse polarization. 

For practical purposes we therefore made the crude but common as- 
sumption, that the transverse asymmetry of the QCD-annihilation graph 
equals the Drell-Yan one. 
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3 Sphinx tt 



In order to estimate the transverse double spin asymmetry in the framework 
of our model we developped a new MC program for transversely polarized 
nucleon-nucleon scattering called Sphinx tt. It is based on the Lund Monte 



Carlo Pythia 5.6 |L9[ and is similar to Sphinx, a MC code for longitudinal 
polarized scattering. 

Important features of Sphinx tt are: 

• Polarization is described correctly in the initial state up to the hard 
partonic reaction. The spin orientations are not followed through the 
final states. String breaking and hadronization in general are assumed 
not to depend on the spin orientation. 

• Two sets of polarized parton distributions are available by now. The 
parametrizations of Aq(x,Q 2 ) by Gehrmann/Stirling and Bartel- 
ski/Tatur Jl3| are taken as transversity distributions at Ql = 4GeV 2 



and are evoluted according to the GLAP-equation (^). The implemen- 
tation of other model parametrizations is easy. 

Transverse polarization is correctly taken care of in the initial state 
showering routine. 

The polarized partonic cross sections for high-p^ jet production, 
prompt photon production and lepton pair production in the Drell- 
Yan model are available. Z° physics is not implemented yet. 

The majority of options of Pythia 5.6 is also available for Sphinx tt. 



4 Results 

Here we present the results of our MC calculations of A TT for the PHENIX 
detector at RHIC. 

Our results are presented for ^Js = 50 GeV and an integrated luminos- 
ity of £ = 200pb _1 as well as = 200 GeV and £ = 800pb _1 . This is 
equivalent to 230 days of beam time with 50% efficiency The assumed inte- 
grated luminosity is larger by a factor of 2.5 compared to the standard value 
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as quoted in |3|, || but is necessary in order to detect an adequate number of 
events ||. 

Taking into account this larger integrated luminosity as well as the 
expected Drell-Yan event rates for the STAR detector |4]] the situation for 
this second detector is the same as presented below. 

We assume that both endcaps will have a polar opening angle of 10° — 
35° with respect to the beam axis and will be able to measure myons with 
1.1 < ?/' < 2.5. The central muon detectors are not very well suited for the 
measurement of A TT since they do not cover the full azimuth. 

MC calculations || indicate competition from charm decay for 
M^- < 4 GeV at yfi = 50 GeV and M^- < 5 GeV at ^ = 200 GeV 
so that we set the minimum lepton pair mass accordingly. We did not set a 
cut to exclude the T production region at 8 GeV < M M + M - < 10 GeV because 
the cross section decreases exponentially with increasing M M + At - anyway. 

Figure [3] shows the differential cross section da/dq T for both cm ener- 
gies. The Drell-Yan graph obviously dominates lepton pair production for 
low qx- The QCD Compton graph is the least important reaction in this 
region. In the high qx region we have the opposite situation. To enhance 
the Drell-Yan signal we therefore required qr < 2 GeV. This leads to an 
improvement of the i^ Dy -factor by AK DY = —0.1 for y/s = 50 GeV and 
AK DY = —0.3 for y/s = 200 GeV in our model. We multiplied the usual 
0(a s )-matrix elements of the hard partonic scattering with the corresponding 
form factors (ITU). Afterwards we fitted the width of the intrinsic momen- 



tum to the data at y/s = 62 GeV taken from |20|. Unfortunately, we got 
(kr) = 800 MeV which is very high in comparison to DIS experiments. On 
the other hand we checked that the asymmetry does not depend on (kr) for 
values smaller than 1 GeV. Our results are rather independent of the details. 
One should also note, that we switched off the general initial state show- 
ering routine because it is only suited for the large qx region and to avoid 
double counting. Bremsstrahlung effects are only taken into account by the 
QCD- annihilation and Compton process. 

Figures |] and [5| show the distribution of events in the plane of the 
Bj0rken variables of the partons which participate in the hard interaction. 
Generally, the x- variables at y/s = 200 GeV have smaller values as the ones at 
y/s = 50 GeV. In both cases, there is one region in the x-x-pl&ne where most 
of the Drell-Yan type events are situated. These regions can be characterized 
by 0.3 <x q < 0.7 and 0.01 < x q < 0.04 at yfs = 50 GeV, and 0.1 < x q < 0.4 
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and 0.002 < x- q < 0.02 at y/s = 200 GeV respectively. 

Fortunately, these ranges are also the ones, where the largest asymme- 
tries are expected since the polaized parton distributions multiplied with x 
have a maximum there. Therefore, it would be interesting to restrict the 
detected events to these regions. Furthermore, the exchange terms in equa- 
tion (p|) would disappear. In order to achieve this one has to calculate xa/b 
from the lepton momenta. Then one must decide which one of the incoming 
partons was the quark and which one was the antiquark. Since the majority 
of events fulfill x q > x q this criterium can be used very successfully for this 
purpose. In only 2%-10% of all cases the x-values are not assigned correctly. 

Of course, it is not possible to reconstruct the x- values of the incoming 
partons of the QCD processes because the outgoing quark or gluon jet is 
not detected. Since the jet momentum is quite small in general we don't 
make a big error by using the Drell-Yan formulas for the reconstruction of 
the Bj0rken-x also for the QCD processes. 

Table [l| shows the event rates for the different cases. At y/s = 200 GeV 
they are about ten to twenty times higher than for ^/s = 50 GeV. This is 
mainly due to the fact that the unpolarized parton distributions rise steeply 
for smaller x. However, according to the shape of the transversity distribu- 
tions for y/s = 200 GeV one should expect a smaller A TT . So we should try to 
minimize the statistical errors by integrating over all kinematical variables. 



The integration over the azimuth is performed as follows [1 I 



D7T 1 — 1 

4 4 '4 

TT ) I f i f f I j / aTT 



■\" { i J ■ J J \*4>A ii y>) ■ (12) 

4 4 — 4 

For cos 20 this gives simply a factor of 4. According to our model the QCD 
annihilation graph carries the same asymmetry as the Drell-Yan graph which 
leads to 

j^TT _ jVpY + jVa nnih , ^TT /-^ 

NbY + -Vannih. + ^Compton 



The statistical error of the asymmetry is 1/yN where iV denotes the 
total event rate. 

Table § contains the results for the asymmetry together with the abso- 
lute statistical error. We can summarize them as follows: 
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If Sq{x,4GeV 2 ) = Aq(x,4GeV 2 ) then A TT won't be measurable at 
RHIC at t/s = 50 GeV because the statistical error is approximately 
as large as the asymmetry itself. 

At t/s = 200 GeV the relative statistical error is also very large, namely 
60% — 80%. Therefore, it is very unlikely that the asymmetry can be 
measured as there will certainly be some unavoidable background and 
systematic error. 

The main reason for the predicted smallness of the asymmetry is the 
very small antiquark polarization which is of the order of 5% accord- 
ing to common parametrizations of parton distributions. By assum- 



ing a much larger sea polarization one gets larger asymmetries [fLC 
Consequently, there is a good chance that A TT can be measured at 
= 200 GeV if, e.g., 5q(x, 4 GeV 2 ) = 3Aq(x, 4 GeV 2 ). 

It is crucial for the experiment that RHIC will be able run with the 
design luminosity in order to have enough events. 



When we were writing this article we received |2T] which describes a 
confinement model calculation of 5q(x,Q 2 ). The authors also calculate A TT 
and their results are of the same order as the results presented in this paper. 
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PHENIX 


subprocess 


^ = 50 GeV 


y/s = 200 GeV 


gg -> fj+fj,- 


5236 


70390 


qq -> 


2157 


34720 




637 


10770 


sum iV 


8030 


115880 


X-f actor 


1.5 


1.6 



PHENIX, x-x-Cuts 



subprocess 



yfe = 50 GeV 



y/s = 200 GeV 



qq -> 
gg -> 
g# -> 



g 



3302 
1027 
353 



58020 
27540 
8740 



sum N 
^-factor 



4682 
1.4 



94300 
1.6 



Table 1: Event rates. 



PHENIX 


parton distribution 


yfe = 50 GeV 


y/s = 200 GeV 


Germann / Stirling 
Bartelski/Tatur 


1.00% ± 1.1% 
0.63% ± 1.1% 


0.38% ± 0.30% 
0.51% ± 0.30% 


PHENIX, x-x-Cuts 


parton distribution 


yfe = 50 GeV 


y/s = 200 GeV 


Germann / Stirling 
Bartelski/Tatur 


1.05% ±1.4% 
1.70% ± 1.4% 


0.45% ± 0.35% 
0.54% ± 0.35% 



Table 2: Transverse double spin asymmetry A TT with absolute statistical 
error. The forseen polarization of beam and target at RHIC (70%) is already 
taken into account. 
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Figure 3: Distribution of the the transverse momentum of the virtual photon 
qr at y/s = 50 GeV and y/s = 200 GeV. qq denotes the Drell-Yan graph, qq B 
the QCD annihilation graph (B stands for Bremsstahlung). 
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Figure 4: Event distribution in the plane of the Bj0rken-a; of the partons that 
participate in the hard reaction at ^/s = 50 GeV. qq denotes the Drell-Yan 
graph, qq B the QCD annihilation graph. 
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Figure 5: Same as figure § for y/s = 200 GeV. 
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